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ABSTRACT 

The  linear  boundary  layer  solution  (Chang,  1973b), 
which  includes  the  temporal  acceleration,  is  utilized 
for  a  CISK  parameterization  to  study  tropical  lower  tro- 
pospheric  waves.   A  one-level  equation  is  derived  for 
the  top  of  the  boundary  layer  and  is  solved  numerically 
for  its  eigenvalues.   Calculations  are  made  for  varying 
boundary-layer  depths  and  vertical  heating  gradients. 
Growth  rates,  periods  and  eigenf unctions  for  the  lowest 
meridional  modes  of  synoptic-scale  disturbances  are 
analyzed.   Results  indicate  that  CISK  heating  does  not 
prefer  any  finite  wavelength. 
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I.    INTRODUCTION 


Observational  evidence  of  tropical  lower  tropospheric 
waves  has  been  extensively  collected   and  studied  during 
the  past  several  years.   Although  there  is  still  some  doubt 
as  to  the  detailed  characteristics  of  these  disturbances, 
it  is  generally  observed  that  these  waves  propagate  west- 
ward, exhibit  a  period  of  about  four  to  five  days  in  their 
meridional  wind  component,  and  may  fall  in  one  of  the  two 
wavelength  ranges:   around  2000-4000  km,  or  approximately 
8000  km.   As  a  direct  result  of  the  observations,  meteoro- 
logists have  attempted  to  explain  theoretically  the 
mechanism  and  structure  of  these  tropical  disturbances. 

Three  different  hypotheses  have  thus  far  been  postu- 
lated:  1)  lateral  forcing  by  mid-latitude  systems,  2) 
barotropic  instability,  and  3)  thermal  forcing.   A  concise 
but  excellent  description  of  these  theories  is  found  in 
Yanai  (1971).   To  summarize  his  discussion  briefly,  various 
investigators  have  found  that  all  three  hypotheses  may  be 
responsible  for  tropical  wave  growth. 

Widespread  attention  has  been  given  the  thermal  forcing 
mechanism  during  the  past  decade.   Ever  since  the  inception 
of  the  so-called  CISK  (Conditional  Instability  of  the 
Second  Kind)  theory,  first  proposed  by  Charney  and  Eliassen 
(1964)  and  Ooyama  (1964)  to  explain  hurricane  growth,  many 


investigators  have  sought  to  extend  its  applicability  to 
the  larger  scale  tropical  waves.   Basically  their  CISK 
models  incorporate  thermal  forcing  by  parameterizing  the 
latent  heat  release  from  small  scale  convection  in  terms 
of  the  boundary  layer  convergence  of  the  large  scale  flows. 
In  addition,  Charney  and  Eliassen's  geos trophic/Ekman 
boundary  layer  formulation  (geostrophic  vorticity  and  ver- 
tical velocity  are  proportional  at  the  top  of  the  boundary 
layer)  was  commonly  used. 

The  primary  approach  in  many  of  the  previous  studies 
was  to  perform  linear  stability  analysis  of  a  system  of 
equations  in  order  to  determine  a  preferred  wave  scale  - 
that  scale  which  has  the  largest  growth  rate  when  other 
conditions  remain  constant.   Most  of  these  earlier  linear 
studies  used  either  a  geos trophic/Ekman  boundary  layer  or 
a  surface  drag  representation  for  the  lowest  level  in  the 
model.  '  Although  the  linearization  process  was  done  to 
simplify  an  analysis,  a  question  generally  arose  in  the 
Ekman-type  boundary  layer  models  near  the  equator.   In  this 
region  the  Coriolis  force  in  the  equations  of  motion 
approaches  zero  and  therefore  the  geostrophic  relation  may 
no  longer  be  valid.   It  is  pointed  out  by  Holton  et  al 
(1971)  that  a  transition  zone  exists  near  the  equator  where 
the  Ekman  boundary  layer  structure  of  the  more  poleward 
latitudes  changes  rather  drastically  into  a  Stokes-type 
boundary  layer  in  the  equatorial  region.   Generally,  this 


transition  zone  has  not  been  accounted  for  by  most  of  the 
investigators . 

For  some  of  the  earlier  CISK  wave  models,  a  major  dis- 
crepancy exists.   Yamasaki  (1971)  and  Murakami  (1972) 
found  that  preferred  wave  scales  on  the  order  of  several 
thousand  kilometers  exist  in  their  models.   However,  in 
contrast  to  these  results,  Chang  (1971)  and  Williams  and 
Robertson  (1973)  showed  that  a  distinct  maximum  growth  rate 
for  the  observed  tropical  waves  did  not  occur. 

Because  of  this  discrepancy,  it  remained  inconclusive 
as  to  the  effect  CISK-type  heating  had  upon  tropical  waves 
in  the  lower  troposphere.   The  purpose  of  this  research  is 
to  use  an  improved  boundary  layer  representation  within  the 
scope  of  linear  theory  to  study  the  effect  of  CISK  upon 
tropical  wave  growth.   The  linear  boundary  layer  solution 
(Chang,  1973b)  which  includes  the  temporal  acceleration  is 
used  for  the  CISK  parameterization.   By  neglecting  the  small 
temperature  fluctuations  in  the  thermodynamic  equation, 
which  is  generally  valid  for  weak  synoptic  -  scale  disturbances 
in  the  lower  tropical  atmosphere,  a  one-level  equation  is 
derived  for  the  top  of  the  boundary  layer.   This  equation  is 
similar  to  the  free  atmosphere  equation  in  Chang  (1973a) . 
An  eigenvalue  problem  is  then  formulated  and  the  growth  rates 
of  the  lowest  meridional  modes  are  obtained  numerically. 


II .    THE  MODEL 

A.    BASIC  EQUATIONS 

The  basic  equations  of  the  model  are  the  linearized 
Boussinesq  primitive  equations  for  an  equatorial  beta- 
plane  : 


3u* 
3t* 


3v* 

3t* 


|£1  +  By.T.  +  K.iV 

3z* 


dp*         0  32v* 


(1) 


(2) 


3u* 
3x* 


3v*    3w* 
9  y  *     3  z  * 


(3) 


where 

x*,y*,z*,t*  =  zonal,  meridional,  vertical  and  time 
coordinates,  respectively 

u*  (x* , y* , z* , t* ) =  zonal  velocity 

v*  (x*  ,y *  ,  z* , t*) =  meridional  velocity 

w*  (x* , y * , z* , t * ) =  vertical  velocity 

p*  (x* , y * , t * ) =  pressure  divided  by  density 

K*  =  vertical  eddy  viscosity  coefficient, 
assumed  constant 

3   =  local  change  of  Coriolis  parameter,  F, 
with  respect  to  latitude,  y,  at  the 
equator;  hence,  equal  to  2ft/a  with  ft 
being  the  angular  velocity  of  the  earth 
and  a  being  the  radius  of  the  earth. 

In  addition,  the  thermodynamic  energy  equation 


3T*      3t*      3t* 


(4) 
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where 


T* 

S* 


temperature 

q      9T* 
static  stability  parameter  {-    - —  + 


where  g  is  the  acceleration 


9z*  ' 


due  to  gravity  and  c   is  the 

P 
specific  heat  of  dry  air  at  constant 

pressure ) 
Q*  =  diabatic  heating  rate 
is  used . 

The  local  change  and  horizontal  advection  terms  in  (4) 
can  be  neglected  considering  the  magnitude  of  each  term 
(Wallace,  1971;  Holton,  1972)  .   The  thermodynamic  equation 
then  becomes 


w*S* 


(5) 


This  equation  implies  that  the  diabatic  heating  is  balanced 
by  the  adiabatic  vertical  motion  field. 

The  basic  equations  are  non-d imens ionalized  using  a 

1/2 
gravity  wave  speed,  c  =  (gH)     ,  where  H  is  the  scale  height 

and  the  depth  scale,  and  the  Rossby  parameter,  3-   These  two 

parameters  are  used  to  derive  the  non-dimensional  units: 

1)  time,  T  =  (c3)~    ,  2)  horizontal  length,  L  =  (c/6)     , 

2 
3)  horizontal  velocity,  c  =  L/T  ,  4)  pressure,  IT  =  c  /L  , 


5)  temperature, 


<J>  =  c  c  /RL 


c  /R,  and  6)  diabatic  heating  rate, 


The  non-dimensional  equations  are 


8u 


3p 


dt       9x 


9^  ^\ 
u 

8z2 


(6) 


11 


3t 


_9p 

9y 


yu  +  K 


9^  /\ 
V 

dz 


u     dv      dw 
dx     dy      dz 


w  S  =  Q  , 


(7) 

(8) 
(9) 


where  the  non-dimensional  static  stability  parameter  is 


RT  dT 

s  =  +  _ 

c  dz 

P 


Solutions  for  (6)  -  (9)  are  assumed  to  be  propagating  waves 
with  zonal  wave  number   A,  and  frequency   V: 


u 

V 

/\ 

w 

p 
Q  ' 


f   u (y , z  )  ^ 
v (y , z) 
w (y , z  ) 

p(y) 

Q(y,z) 


ei(Ax+Vt,.(10) 


ivu  =  -iXp  + 


32 

d  U 
p  +  yv  +  K  -  , 


Substituting  (10)  into  (6)  -  (9),  the  equations  become 

(11) 
(12) 
(13) 


1W  =  - 


3y 


s2 

yu  +  K   —  , 


8z 


.  1       dv     8w 
iAu  +  t—  +  x—   =  0  , 
dy    dz 


W  S  =  Q 


(14) 
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By  eliminating  u,  v  and  w  in  (11)  -  (14) ,  a  single 
differential  equation  in  p  results: 


3y2 


"       [ 


2y 


2-2 

y       v 


9p 

Sy 


+     [ 


X    (y    +   v    ) 


(y2- 


v2, 


A2]    p 


My 


-   v2, 


V 


9z     V 


(15) 


In  this  equation  the  only  vertically  differentiated  term 
is  the  heating  term  on  the  right-hand  side.   This  implies 
that  the  system  of  equations  is  vertically  coupled  only 
through  diabatic  heating.   Hence,  if  the  vertical  profile 
of  Q  is  specified,   (15)  can  be  applied  as  an  ordinary  dif- 
ferential equation  with  y  as  the  only  independent  variable 
at  any  level  above  the  boundary  layer. 


B.    THE  CISK  EQUATION 

The  CISK  process  is  parameterized  so  that  the  diabatic 
heating  function,  Q(z),  is  related  to  the  vertical  velocity, 
w  ,  at  the  top  of  the  boundary  layer  (z  ) : 


Q (z)  -  n (z  )  *  w, 


(16) 


where  n,  ( z )  is  a  proportionality  function  which  specifies 
the  vertical  profile  of  heating.   Notice  that  (16)  also 
implies  negative  heating  when  downward  motion  occurs  at 
the  top  of  the  boundary  layer.   This  is  difficult  to  avoid 
due  to  the  separation  of  the  x-dependence  of  the  solutions 

The  vertical  motion  at  the  top  of  the  boundary  layer 
is  expressed  by 
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S  p       9p    , 

wt  =  Ri  77  +  R2  af  +  AR3p 

3y 


(17) 


where 


R,  =  B  C   -  B  C 

1     +  +     -  - 


R 


2    2KY    +  + 


(B  D   +  B  D  ) 


+  B  C  [(y+V)  *  -  2iy2m  B  A  2] 


-1  2 

+  B_C_[(y-V)     +  2iaz   B_A   ] 


R 


3     2KY    +  + 


(B  D   -  B  D  ) 


and 


+  B+C+[(y+V)~1  -  2iYzT  B+A+2+  A] 


+    B_C_[(y-V)  -     2iaz        B_A     2-     A] 


exp  (-YzT) 


A 
B 
B 
C 
C 
D 


exp ( -az    ) 

|(y+v)_1     (1    -    A+2)_1 

i(y-v)"1     (1    -    A_2)_1 

zm    -    —    +    -    A        (z       +    — 

T  Y  Y  +  T  Y 


1 
z       -    —   + 
T        a 


2A, 

a 


-     A     2(z        +     i) 

-         T        a 


-       2A    -1-A        )     +     2z     A       -    1-A        -     z     A 

2  +  +  T    +    Y  +  T    + 

Y  ' 


12  1 

D       =    -     (2A    -1-A        )  +     2z    A     [-(1-A     )-     z     A     ] 
—               2            ~               ~  T     —     CX  —  T     — 

a 
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and 


2  =  i  (V  +  y) 
K 


a 


2    i(v-y) 
K 


Equation  (17)  is  derived  by  integrating  the  total  boundary 
layer  convergence,  using  the  solutions  in  Chang  (1973b). 
This  boundary  layer  solution  includes  the  temporal  accelera- 
tion and,  therefore,  is  frequency  dependent.   If  the  bound- 
ary layer  depth  becomes  infinite  (z   •>  °°)  ,     a  singularity  in 
w   exists  at  the  critical  latitude  where  y  =  V  (Holton  et  al, 
1971).   If  z       is  finite,  the  solution  may  possess  quasi- 
singular  behavior  around  this  latitude  depending  upon  z 
and  the  mode  of  symmetry  about  the  equator  (Chang,  1973b). 

The  heating  profile,  r\  ( z )  ,  is  specified  such  that  8r|/3z 
varies  from  IS  to  2S  at  the  top  of  the  boundary  layer.   Con- 
sequently (15)  can  then  be  applied  at  z  .    [Notice  that 
(13)  and  (16)  require  that  n,  =  S  at  the  top  of  the  boundary 
layer.   For  a  lapse  rate  of  6°C  km    and  w  =1  cm  sec    the 
magnitude  of  T)    in  turn  produces  a  heating  rate  of  approxi- 
mately 3°C  day   .]   Because  the  slope  of  the  heating  func- 
tion is  defined,  it  is  implied  that  the  heating  function 
varies  smoothly  in  the  vicinity  of  z  =  z  .   Furthermore,  a 
positive  slope  implies  a  heat  source  above  the  top  of  the 
boundary  layer.   Finally,  to  reconcile  the  previous 


15 


assumption  that  the  pressure  perturbation  is  independent 

of  height  for  z  <  z  ,  it  is  assumed  that  the  smaller  but 

—   T 

finite  heating  in  the  uppermost  layer  of  the  boundary 
layer  acts  in  such  a  way  that  p  remains  constant  in  height 

A  second-order  homogeneous  differential  equation  in  p 
is  obtained  by  substituting  (16)  and  (17)  into  (15) : 


7*  +p   If +  =  p  =  °  - 

dy 


(18) 


where  E,  F,  and  G  are  functions  of  y  and  v 


E     =     1     + 


F     =    - 


■  ,     2  2, 

i  (y    -    v    ) 


vs 


R. 


dz 


2y 


2  2 

(y   -    v    ) 


1  (y    -    v    ) 

vs 


R 


dr, 
dz     ' 


_   _    A         (y2+   v2)         ,2  i(y2-   v2)         Xt?  dn 

G  =  v     ~~2     2~  '  A    +         ^rr       AR3     dz 

(y    -   v    ) 
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Ill .    SOLUTION  PROCEDURE 

To  solve  (18)  numerically,  a  central  difference  scheme 
in  terms  of  the  independent  variable,  y,  is  utilized.  The 
difference  expression  is 


(Ei+Fi  !*>  pi+i+[Gi  (Ay)2-2  Ei]  v(Ei-Fi  f*>  pi-i =  °  (i9) 


i  =  2,  3,  4,...,  J-l, 


where  J  is  the  total  number  of  grid  points  and 


Ay  = 


width  between  the  two  boundaries 
J-l 


A  channel-flow  boundary  condition  is  assumed  such  that  the 
perturbations  vanish  at  the  boundaries;  that  is, 


Pj  =  0 


(20) 


The  difference  equation,   (19) ,  is  solved  for  its  eigen- 
value, V,  using  an  iterative  procedure  as  described  by 
Stone  (1970)  .   First,  for  a  given  value  of  X,  z  ,     K  and 
8r|/8z  ,  an  initial  guess  for  v  is  made.   Then  by  satisfying 
one  of  the  boundary  conditions  in  (20) ,  the  pressure  per- 
turbation at  the  first  grid  point  is  set  to  zero;  that  is, 
p   =  0.   Next,  the  amplitude  of  the  eigenf unct i on  at  the 
second  grid  point  is  arbitrarily  chosen  as  unity.   Using 
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these  first  two  values  for  the  e igenf unctions ,  the  numeri- 
cal solution  process  successively  solves  for  the  remaining 
values  of  the  eigenf unctions :   p  ,  p.,...p  . 

In  general,  the  guess  for  v  is  incorrect  and  the  calcu- 
lated value  for  p   is  therefore  non-zero,  thereby  violating 

J 

the  second  boundary  condition  of  (20) .   Consequently,  the 
calculation  of  p   is  repeated,  using  a  value  of  v  which  dif- 
fers  from  the  initial  guess  by  1%.   The  two  values  of  p   are 
then  used  to  perform  a  linear  extrapolation  from  the  two 
valines  of  V  to  a  new  value  of  V  which  would  make  p   zero  if 
it  depends  on  V  in  a  linear  fashion.   The  new  value  of  V  is 
then  used  to  reiterate  the  whole  procedure,  and  iterations 
are  repeated  until  a  value  of  v  accurate  to  0.01%  is  ob- 
tained.  In  summary  the  process  is  as  follows: 


N+l     N-l 
V     =  V 


Av 
Ap 


N-l 


(21) 


where  N  is  the  iteration  step, 


A       N     N-l 

Av  =  V   -  V 


.         N      N-l 
APJ  =  PJ  "  PJ 


The  upper  limit  for  N  is  prespecified  as  25  in  almost  all 
the  computer  runs.   Convergence  to  a  solution  generally 
occurs  well  before  this  upper  bound  but,  if  convergence  is 
not  achieved  by  N=25,  that  particular  run  is  terminated 
and  a  new  set  of  input  values  is  automatically  provided  to 
maximize  computer  efficiency. 


18 


z       =     1    km    to    4    km    and    8f|/9z    =     IS    to    2S.        The    solu> 


Two  external  parameters  are  varied  in  the  following 
range : 

tions  are  obtained  for  both  the  symmetric  (pressure  is 
symmetric  about  the  equator)  and  the  asymmetric  (pressure 
is  asymmetric  about  the  equator)  cases. 
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IV.    RESULTS 

A.    FREE  MODE 

Solutions  for  the  "free"  mode,  when  the  right-hand 
side  of  (i5)  is  set  to  zero,  are  obtained  as  a  test  for 
the  model.   Zero  growth  rates  (infinite  e-folding  times) 
result  for  all  wavelengths. 

In  Figure  1,  the  results  for  wave  period  are  shown. 
Notice  that  wavelengths  of  the  synoptic  scale  range  (2000 
-10,000  km)  have  periods  between  approximately  three  and 
seven  days  for  the  asymmetric  case  and  five  and  ten  days 
for  the  symmetric  case.   In  general,  the  longer  waves 
have  higher  frequencies  than  the  shorter  waves  in  both 
cases.   These  free  mode  frequencies  correspond  quite 
closely  to  those  described  by  the  barotropic  Rossby  wave 
frequency  formula 


V  =  A/(A2+  I2)     , 


where  i    =  tt/D,  D  being  the  width  of  the  channel,  and  the 
approximated  Matsuno  (1966)  wave- frequency  formula 

V  =  k/(k2+  2n  +  1)   . 

Here  k  is  the  non-dimensional  zonal  wave  number  and  n  is 
the  meridional  mode  number. 
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B.    GROWTH  RATE 

Figures  2  through  5  indicate  the  growth  rates  for  both 
symmetry  cases.   The  important  result  shown  by  all  four  of 
these  figures  is  that  the  synoptic  scale  waves  have  essen- 
tially the  same  growth  rate.   Thus,  there  exists  no  pre- 
ferred wavelength  in  the  synoptic  scale  range  -  a  result 
which  agrees  with  those  of  Chang  (1971)  and  Robertson  (1973) 

For  each  set  of  external  parameters,  two  different  solu- 
tions are  generally  found:   a  long-wave  (LW)  mode,  one  that 
possesses  greater  growth  rates  in  the  long-wave  portion  of 
the  synoptic  scale;  and  a  short-wave  (SW)  mode,  one  that 
favors  short-wave  growth.   Both  modes  of  the  symmetric  case 
generally  show  larger  growth  rates  (faster  e-folding  times) 
than  those  of  the  asymmetric  case,  although  for  shorter 
wavelengths  this  difference  rapidly  decreases. 

An  increase  in  the  slope  of  the  heating  function,  r\  , 
increases  the  growth  rates  of  the  waves.   For  a  boundary 
layer  depth  of  one  kilometer,  both  modes  of  the  asymmetric 
and  symmetric  cases  indicate  an  approximate  doubling  of 
their  growth  rates  with  an  increase  of  Br)/9z  from  Is  to  2s 
(compare  Figures  2  and  4) .   For  z   =2  km,  the  same  increase 
in  the  vertical  gradient  of  heating  produces  an  even  larger 
change  in  the  growth  rates  -  generally  a  125%  increase 
results  (see  Figures  3  and  5) . 

Likewise,  a  deeper  boundary  layer  produces  larger  growth 
rates.   This  effect  is  especially  noticeable  for  the  LW 
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mode  of  the  asymmetric  case.   Growth  rates  for  this  mode 
approximately  triple  when  the  boundary  layer  deepens  from 
one  to  two  kilometers,  but  increase  more  slowly  with 
deeper  boundary  layers  (results  not  shown) .   The  SW  mode 
for  the  asymmetric  case,  as  well  as  the  two  symmetric 
case  modes,  exhibits  similar  but  less  pronounced  increases 
in  growth  rates  when  the  boundary -layer  depth  increases 
from  one  to  two  kilometers. 

A  rather  interesting  result  is  that  the  symmetric  SW 
mode  for  shallow  boundary -lay  er  depths  (z   <  1.3  km  for  a 
wavelength  of  4500  km) ,  a  low  vertical  heating  gradient 
(dr\/d?-    -    IS)  and  longer  (^4000  km)  wavelengths  has  larger 
growth  rates  than  the  LW  mode  (Figure  4) .   But  for  deeper 
boundary  layers  the  opposite  occurs;  that  is,  the  growth 
rate  of  the  symmetric  LW  mode  is  greater  than  those  of  the 
SW  mode  (Figure  5) .   Thus,  a  simple  relationship  between 
the  growth  rates  of  the  symmetric  modes  and  the  boundary- 
layer  depth  is  not  apparent.   On  the  other  hand,  Figures 
2  and  3  indicate  the  asymmetric  SW  mode  has  larger  growth 
rates  than  those  of  the  LW  mode  for  both  bounda ry- layer 
depths . 

C.   PERIOD 

As  was  mentioned  earlier  the  periods  of  the  free  modes 
(Figure  1)  were  determined  as  a  test  case.   Figures  6 
through  9  illustrate  the  effects  of  heating  and  boundary- 
layer  depth  upon  wave  period.   Note  that  the  periods  for 
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the  LW  modes  of  both  the  asymmetric  case  and  symmetric 
cases  very  closely  resemble  those  of  the  corresponding 
free  modes.   This  result  is  expected,  because  the  LW  modes 
generally  have  much  smaller  growth  rates  than  the  SW  modes, 
thereby  better  approximating  the  free  modes'  zero  growth 
rates.   Both  the  LW  and  SW  modes  of  the  asymmetric  case 
exhibit  shorter  periods  than  those  of  the  corresponding 
modes  of  the  symmetric  case. 

It  is  also  noted  that  varying  3  H  /  9  z  produces  little 
change  in  frequency.   Only  a  slight  decrease  in  period  is 
noticeable  with  an  increase  in  9r|/9z.   The  deeper  boundary 
layer  does  not  appreciably  affect  the  wave  period  either, 
although  a  very  small  increase  in  frequency  is  detectable 
for  the  longer  wavelengths  of  the  symmetric  LW  mode  and 
a  small  decrease  for  the  asymmetric  LW  and  SW  modes,  as 
well  as  for  the  symmetric  SW  mode. 

D.    EIGENFUNCTIONS 

Figures  10  -and  11,  which  represent  the  free  mode  pro- 
files for  z   =1  and  2  km,  respectively,  are  consistent 
with  the  solutions  in  Chang  (1973b) .   When  CISK  heating 
is  included,  the  p  and  w   profiles  for  both  symmetry  cases 
generally  tilt  northeast-southwest  between  the  northern 
boundary  and  the  equator.   However,  the  equatorward  limit 
for  the  tilt  in  the  p  profile  of  the  symmetric  case  depends 
upon  wavelength  (the  longer  the  wavelength  the  farther  away 
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from  the  equator  is  the  limit).   Such  a  northeast-southwest 
tilt  in  the  eigenf unctions  implies  poleward  transport  of 
westerly  momentum  and  agrees  with  observations  summarized 
in  Wallace  (1971) . 

Comparison  of  Figures  12  and  16  (among  others)  shows 
that  the  tilt  for  the  asymmetric  LW  mode  decreases  with 
longer  waves.   The  SW  mode  generally  shows  a  similar  change. 

By  increasing  the  vertical  heating  gradient  at  the  top 
of  the  boundary  layer,  the  tilt  increases  for  both  asymme- 
tric modes  especially  in  the  short  waves  (see  Figures  12- 
15) .   The  only  significant  change  in  tilt  due  to  an  increase 
in  the  boundary- lay er  depth  occurs  for  the  LW  mode,  for 
which  the  tilt  increases  with  decreasing  wavelength  (compare 
Figures  12  with  14  and  13  with  15). 

Figures  12  and  16  show  a  decrease  in  the  SW  mode's 

equatorward  maximum  of  w   amplitude  and  a  general  smoothing 

of  the  w   profile  when  the  wavelength  decreases.   On  the 

other  hand,  the  LW  mode  shows  a  slight  equatorward  shift 

of  the  w   maximum  and  a  uniform  amplitude  increase  with 
T  ^ 

decreasing  wavelength.   The  p  profiles  for  both  modes  re- 
main nearly  unchanged  in  these  figures. 

An  increase  in  df\/dz    shifts  the  asymmetric  LW  mode's 

maximum  in  w   and  p  more  equatorward,  particularly  for  the 
T 

shorter  waves  (see  Figures  12  and  13) .   The  SW  mode,  mean- 
while, shows  a  slight  equatorward  shift  only  in  its  pole- 
ward w   maximum  and  a  somewhat  smoother  profile  of  w   when 
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9r)/3z  increases.   The  p  profile  of  the  SW  mode  remains  al- 
most unchanged  with  a  larger  8r|/9z. 

Deepening  the  boundary  layer  also  effects  an  equator- 
ward  shift  in  the  asymmetric  LW  mode's  p  and  w   maxima 
(compare  Figures  12  with  14  and  16  with  18) .   This  result 
is  similar  to  that  found  by  Chang  (1973b).   From  Figures 
16  and  18,  it  is  apparent  that  for  longer  waves  the  SW 
mode's  equatorward  maximum  in  w   becomes  the  primary  maxi- 
mum when  the  boundary- lay er  depth  increases.   While  both 
w   maxima  of  the  SW  mode  shift  slightly  toward  the  equator 
with  a  deeper  boundary  layer,  the  p  profile  remains 
virtually  unchanged  (Figures  16  and  18) .   Both  modes  indi- 
cate a  larger  w   amplitude  when  z   increases  from  1  km  to 
2  km,  while  the  p  amplitude  again  remains  nearly  constant 
(compare  Figures  12  with  14  and  16  with  18) . 

Essentially  the  same  results  as  those  just  noted  for 
the  asymmetric  case  are  found  for  the  symmetric  case. 
Only  the  significant  differences  are  described  below. 

From  Figures  24  and  25  it  is  seen  that  the  equatorward 
w   maximum  of  the  symmetric  SW  mode  shifts  poleward  for 
the  longer  wavelengths  when  9n,/9z  increases.   In  contrast, 
the  asymmetric  SW  mode's  equatorward  w   maximum  remains 
stationary  when  9ri/9z  is  increased  but  the  poleward  maxi- 
mum moves  equatorward.   While  an  increase  in  the  boundary- 
layer  depth  shifts  the  asymmetric  LW  mode's  w   and  p  maxima 

more  toward  the  equator,  only  the  w   maximum  moves 

T 
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equatorward  for  the  symmetric  LW  mode  (compare  Figure  20 
with  22  and  24  with  26). 

The  phase  difference  between  the  pressure  and  vertical 
motion  fields  generally  remains  1/2  cycle  for  both  symme- 
try cases  under  the  conditions  investigated.   The  only 
exception  appears  with  the  symmetric  case.   In  the  vicinity 
of  the  critical  latitude  and  equatorward,  a  rapid  phase 

shift  in  wm  occurs  such  that  the  two  fields  are  nearly  in 

T  J 

phase  at  the  equator.   Whereas  an  out-of-phase  relation- 
ship implies  a  downward  energy  transfer  from  the  source 
region  immediately  above  the  boundary  layer,  in-phase  im- 
plies the  converse  transport  of  energy,  that  is,  from  the 
boundary  layer  to  the  atmospheric  layer  immediately  above. 
However,  the  relative  amplitudes  of  both  p  and  w   are 
small  near  the  equator,  thereby  minimizing  the  upward 
transport  of  energy  within  the  near -equator ial  region. 
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V.    SUMMARY  AND  CONCLUSIONS 

A  linear,  tropical  boundary- lay er  model  is  formulated 
on  an  equatorial  beta-plane  by  utilizing  the  primitive 
equations  of  motion,  a  simplified  thermodynamic  energy 
equation  and  a  CISK  parameterization.   Solutions  in  the 
form  of  complex  wave  frequencies  are  obtained  numerically 
as  eigenvalues  to  a  second-order  finite  difference  equa- 
tion . 

The  primary  conclusion  reached  is  that  there  is  no 
preferred  synoptic-scale  wavelength  for  CISK-forced  tropi' 
cal  waves.   Although  two  modes  for  each  symmetry  case 
exist  when  the  CISK  forcing  is  incorporated.,  in  no  case 
does  a  noticeable  peak  in  the  growth  rate  curves 
materialize  in  the  synoptic  scale  range.   In  addition, 
the  e-folding  times  for  the  more  realistic  vertical  heat- 
ing gradient  (dr\/dz    =    Is )  are  found  to  be  of  the  order  of 
5  days  or  slower,  even  for  the  deeper  boundary- lay er 
depths.   This  magnitude  of  growth  rate  is  quite  small 
compared  to  baroclinical ly  unstable  waves  in  the  midlati- 
tudes.   Consequently,  even  if  one  scale  is  slightly  pre- 
ferred among  the  others,  its  relatively  faster  growth 
would  be  at  such  a  slow  rate  that  the  competition  for 
energy  among  the  waves  may  not  be  strong  enough  to  pro- 
duce a  dominant  wavelength. 
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The  main  conclusion  of  this  study  -  that  the  growth 
rate  is  nearly  uniform  in  the  synoptic  scale  range  -  agrees 
with  some  earlier  works  (Chang,  1971, and  Williams  and 
Robertson,  1973),  but  differs  from  others  (Yamasaki,  1969, 
1971;  Hayashi,  1971;and  Murakami,  1972).   A  similar  inde- 
pendence of  scale  is  noted  for  the  so-called  "spin-down" 
problem.   By  specifying  dr\/dz    =  -1  and  substituting  it  into 
the  forcing  function  on  the  right-hand  side  of  115) ,  a 
crude  formulation  of  the  spin-down  problem  is  made.   By  so 
doing,  the  vertical  motion  profile  for  the  free  atmosphere 
is  assumed  to  decrease  linearly  from  w   at  the  top  of  the 
boundary  layer  to  zero  at  one  scale  height  (i.e. ,  z  =  H) . 
The  resultant  damping  rates  are  found  by  using  the  same 
procedure  for  computing  the  CISK  growth  rates  and  they  are 
shown  in  Figure  28.   The  spin-down  times  are  small  as 
expected  because  of  the  small  rotation  rates  near  the 
equator.   In  addition  the  curves  in  Figure  28  show  nearly 
uniform  spin-down  rates  for  the  synoptic-scale  waves.   The 
CISK  model  equation,  therefore,  is  mathematically  analogous 
to  the  spin-down  equation  except  that  the  forcing  function 
has  the  opposite  sign.   Hence,  it  may  not  be  surprising 
that  the  CISK  growth  rates  are  also  independent  of  scale. 

Although  there  is  no  preferred  wave  scale  in  the  2000- 
10,000  km  range,  the  LW  mode  of  both  symmetry  cases  never- 
theless exhibits  larger  growth  rates  in  the  synoptic  scale. 
In  contrast,  the  short  wave  end  of  the  LW  mode  has  decreas- 
ing growth  rates  with  decreasing  wavelength.   This  mode 
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with  its  growth  rate  characteristics  seems  to  agree  with 
the  HB  mode  of  Yamasaki  (1969)  and  the  findings  of  Murakami 
(1972).   In  addition,  their  results  are  similar  to  the  LW 
mode  in  the  following  two  ways:   1)  maximum  amplitude 
exists  near  the  equator,  and  2)  wave  periods  of  about  5 
days  occur.   On  the  other  hand,  the  SW  mode,  especially 
for  the  shallower  boundary  layers,  has  its  maximum  ampli- 
tude more  poleward  and  possesses  larger  periods.   The 
SW  mode  also  appears  to  prefer  zero  scale  for  growth. 
These  three  characteristics  indicate  that  the  SW  mode  may 
be  identified  as  a  hurricane  (or  typhoon)  mode. 

From  the  results  of  this  study,  a  growth  rate  depen- 
dence upon  boundary- lay er  depth  appears  to  be  two-fold. 
First,  by  increasing  the  depth  of  the  boundary  layer, 
growth  rates  increase  for  both  the  LW  and  SW  modes  of  both 
symmetry  cases.   Secondly,  the  greatest  change  in  the 
growth  rates  for  both  cases  generally  occurs  between  z 

=  1  km  and  z   =2  km.   As  a  result  of  these  findings,  it 
T 

is  apparent  that  linear,  tropical  boundary- layer  models 
must  account  for  the  variations  in  boundary- lay er  depth 
in  order  to  describe  more  accurately  the  growth  rates  of 
tropical  lower  tropospheric  waves  under  CISK  forcing. 

Within  the  tropical  boundary  layer,  there  is  generally 
a  transition  zone  between  a  quasi-Ekman  boundary  layer  of 
the  more  poleward  latitudes  and  a  quasi-Stokes  boundary 
layer  of  latitudes  near  the  equator.   The  slow-rotation 
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Stokes  boundary  layer  with  its  inefficient  pumping  does  not 
contribute  much  to  wave  growth.   The  inclusion  of  the  tem- 
poral acceleration  terms  in  the  bou ndary- lay er  equations 
thus  seems  to  have  very  little  effect  in  CISK  instability. 
In  fact  the  results  of  this  study  generally  agree  with  the 
quasi-geostrophic  models  of  Chang  (1971)  and  Williams  and 
Robertson  (1973) . 

At  the  top  of  the  boundary  layer,  the  vertical  motion, 

w  ,  assumes  a  smoother  profile  across  the  transition  zone 
T 

(critical  latitude)  for  the  growing  solution  compared  to 
the  free  solution  in  Chang  (1973b) .   However,  the  large 
phase  shift  across  the  transition  zone  remains  similar  to 
that  found  for  the  free  solution. 

Non-linear  advective  terms  may  be  important  in  affect- 
ing the  growth  rates  of  synoptic- scale  tropical  waves. 
Such  terms  are  neglected  in  this  study  and  should  be  in- 
cluded in  future  research  to  determine  their  effect  in  a 
CISK  model.   Finally,  a  possible  improvement  to  this  linear 
model  may  be  the  inclusion  of  lateral  shear  of  the  mean 
flow  in  the  equations  for  both  the  boundary  layer  and  the 
free  atmosphere. 
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Figure    1.       Wave    periods    for    free    modes. 
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Figure    2.       Growth    rates    for    asymmetric    case    and    z       =     1    km 
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Fiqure  3.   Growth  rates  for  asymmetric  case  and  z   =2  km 
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Figure  4.   Growth  rates  for  symmetric  case  and  z   =  1  km 
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Figure  5.   Growth  rates  for  symmetric  case  and  z T  =  2  km 
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Figure  6.   Wave  periods  for  asymmetric  case  and  z T  -  !  km 
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Figure    7.       Wave    periods    for    asymmetric    case    and     z^    -    2    km 
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Figure 


Wave  periods  for  symmetric  case  and  z   =  1  km 
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Figure  9.   Wave  periods  for  symmetric  case  and  z   -  2  km 
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Figure    11 


Eigenf unctions  for  free  mode  of  symmetric 
case  (two  wavelengths  given) .   Solid  lines 
indicate  pressure  perturbation,  p,  and 
dashed  lines  the  vertical  motion,  w  ,  at  the 
top  of  the  boundary  layer. 
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Figure  12.   Eigenf unctions  for  2000  km  wavelength, 

asymmetric  case,  z  =  1  km  and  8n,/3z  =  IS. 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  13.   Eigenf unctions  for  2000  km  wavelength, 

asymmetric  case,  z=    1  km  and  8  n/  3  z  =  2S  . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  14.   Eigenf unctions  for  2000  km  wavelength, 

asymmetric  case,  z  =  2  km  and  dr)/dz    =    IS. 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  .  at  the  top  of  the  boundary  layer. 
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Figure  15.   Eigenf unctions  for  2000  km  wavelength, 

asymmetric  case,  z  =  2  km  and  8r)/8z  =  2S • 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  16.   Eigenf unctions  for  4500  km  wavelength, 

asymmetric  case,  z  =  1  km  and  dr\/dz    =    IS. 
Solid  lines  indicate  pressure  perturbation, 


T' 


and  dashed  lines  the  vertical  motion, 
it  the  top  of  the  boundary  layer. 
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Figure    17 . 


Eigenf unctions  for  4500  km  wavelength, 
asymmetric  case,  z   =  1  km  and  8t]/8z  =  2S  . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 


w. 


,t  the  top  of  the  boundary  layer 
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Figure  18.   Eigenf unctions  for  4500  km  wavelength, 

asymmetric  case,  z    =    2    km  and  dr)/dz    =    Is  . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  19.   Eigenf  uric  t  ions  for  4500  km  wavelength, 

asymmetric  case,  z    =    2    km  and  dr\/dz    =    2S  . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure    20. 


Eigenfunctions  for  2000  km  wavelength, 

symmetric  case,  z    -    1  km  and  3n/9z  =  IS. 

Solid  lines  indicate  pressure  perturbation, 

p,  and  dashed  lines  the  vertical  motion, 

w  ,  at  the  top  of  the  boundary  layer. 
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Figure  21.   Eigenf unctions  for  2000  km  wavelength, 

symmetric  case,  z  =  1  km  and  8r]/3z  =  2S  . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  22.   Eigenf unctions  for  2000  km  wavelength, 
symmetric  case,  z_=  2  km  and 


Solid  lines  indicate  pressure  perturbation, 
and  dashed  lines  the  vertical  motion, 
at  the  top  of  the  boundary  layer. 


52 


AMPLITUDE 


PHASE 


Figure  23.   Eigenf unc t ions  for  2000  km  wavelength, 

symmetric  case,  z  =  2  km  and  dr]/dz    =    2S  . 
Solid  lines  indicate  pressure  per truba tion , 
p,  and  dashed  lines  the  vertical  motion, 
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T' 


t  the  top  of  the  boundary  layer 
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Figure  24.   Eigenf unct ions  for  4500  km  wavelength, 

=  1  km  and  8n/9z  =  is 


case 


symme tr ic 

Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  25.   Eigenf unctions  for  4500  km  wavelength, 

symmetric  case,  z  =  1  km  and  8n/3z  =  2S . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure  26.   Eigenf unctions  for  4500  km  wavelength, 

symmetric  case,  z  =  2  km  and  8n./3z  =  Is  . 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  ,  at  the  top  of  the  boundary  layer. 
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Figure    27.       Eig 


iqenf unctions  for  4500  km  wavelength, 


symmetric  case,  z  =  2  km  and  8n/8z  =  2S. 
Solid  lines  indicate  pressure  perturbation, 
p,  and  dashed  lines  the  vertical  motion, 
w  .  at  the  top  of  the  boundary  layer. 
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Figure    28.        Spin-down    times    for    z       =    2    km. 
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